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ABSTRACT 

We study the prospects for measuring the dark matter distribution of voids with 
stacked weak lensing. We select voids from a large set of N-body simulations, and 
explore their lensing signals with the full ray-tracing simulations including the effect 
of the large-scale structure along the line-of-sight. The lensing signals are compared 
with simple void model predictions to reconstruct the three-dimensional mass dis- 
tribution of voids. We show that the stacked weak lensing signals are detected at 
significant level (S/N^ 5) for a 5000 degree 2 survey area, for a wide range of void 
radii up to ~ 50 Mpc. The error from the shape noise little affects lensing signals at 
large scale. It is also found that dense ridges around voids have a great impact on 
the weak lensing signals, suggesting that proper modeling of the void density pro- 
file including surrounding ridges is essential for extracting the average total mass of 
voids. 

Key words: gravitational lensing: weak, stacked lensing, large-scale structure of uni- 
verse 
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1 INTRODUCTION 

Cosmological structure formation and matter distribution of 
the universe depend on models of cosmology. Therefore one 
of the keys for constraining cosmology is to understand the 
distribution of matter. High density peaks in the matter den- 
sity distribution which correspond to galaxies and galaxy 
clusters have been well studied dBroadhurst et alJ 120051; 
lUmetsu et ai]|201ll: lOguri et alj|2012h . These o yerdensity re- 
gions are used for constraining cosm ology iTakadal 120061; 
lOguri & Takadall201ll:rPrada et al.ll2012h . On the other hand, 
low density peaks which correspond to voids have attracted 
less attention in spite of its long history ( Kirshner et aljl98ll) 
and the fact that they fill more than 80% of the volume of the 
universe and therefore are the essential component of large 
scale structure (LSS) of the universe. 

There have been severel analytica l works that studied 
the f o rmation and evolution of voids ^Fillmore & Goldreichl 
11984; ISuto et al.1 11984; iFurlanetto & Piranl 120061) . The early 
evolution of void is well described by a spherical collapse 
model. In this model underdense regions expand a nd grad- 
ually become deep until the shell crossing occurs l lPeeblesI 
1980). Effects of different cosmology on voids have also 
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20071; iLee & Parld 
13). These papers 



been studied jMoffatl 120061 iPark & Led 12007 
l2009HKamionkowski et al.ll2009l;lLi et alj|2012 
pointed out possibilities for constraining cosmology with 
voids, including primordial non-gaussianity and modified 
gravity theories. 

On the other hand, observationally various void find- 
ing algorithms have been applied to redshift surveys such 
as the Slo an Digital Sky Survey (SPSS) to study properties 
of voids iStrauss et al.lll992l; llones et alj|2004; IColberg et all 
l2008l:|Pan et alj2012l) . It was found that voids are significantly 
underdense in their observations which use galaxy distribu- 
tion for searching voids. These results are broadly consistent 
with theoretical predictions. 

However many problems still exist when comparing 
simulations with observations (e.g. iRyden & T urne 3 11984; 
iKlypin et a!lll999l;lTavasoli et alj|2012h . In addition, the mat- 
ter distribution of voids is estimated by assuming the galaxy 
bias on observations, and therefore the uncertainty of the 
bias is one of the most serious challenges in void studies. 
By contrast, gravitational lensing trace all matters regardless 
of baryon or dark matter. In fact, the possibility of study- 
ing voids with weak gravit ational lensing was discussed in 
iTain & Van Waerbeke i20 0j^), a nd a possible detection was 
made in Mivazaki et al. 12002 ). While lensing si gnals of in- 
dividual voids are weak l lAmendola et al.l ll999), next gen- 
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eration wide field s urveys such as Sub aru Hyper Suprime- 
cam survey (HSC) ( Mivazaki et al. 2006) will be able to detect 
such signals and study the true matter distribution in voids 
without any assumption of galaxy bias. 

This paper is organized as follows. In Sec|2] we describe 
basics of our analyzing techniques, focusing on how we can 
obtain information of dark matter distribution from gravita- 
tional lensing. In Sec|3j we describe our simulation and void 
finding algorithm. In Sec|U we describe the mass function 
of voids. The void model used for fitting lensing signals are 
presented in Sec|5] In Sec|6] we show the result of the stack- 
ing analysis and fitting with our model. We summarize our 
results in Sec[7l 
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Figure 1. Definition of the coordinate system 



2 BASICS OF ANALYSIS 

2.1 basics of gravitational lensing 

Throughout this work, we assume a spatially flat universe. 
Gravitational lensing effects on the 0L-f> plane is character- 
ized by an isotropic stretching called convergence jc(0) and 
anisotropic distortion called shear 71 (0) and 72(6). They are 
related to the two dimensional analogue of the Newtonian 
gravitational potential ip(d) as 

K (0) = iv 2 ^(0), (1) 

*<•> - 

Then convergence at the position is represented as 

K{e) = 3Wn m0 D h D, r dzmz)) (4) 
^~c LJgaj J 

where z is distance from lens which are connected with dis- 
tance r from the center of a lens and distance Dj 1 on the lens 
plane (figO 

r 2 = (D!0) 2 + z 2 , 

e 2 = e 2 + e 2 p, 

and 3(8, z) is density contrast 

with p(z) being the average density at z. Hq, O m o are nub- 
ble parameter and matter density parameter at present, c is 
speed of light, aj is scale factor at lens position. D s , Dj and 
Dj s are angular diameter distance which are expressed as 

D s = & s xs, Di = &\Xh Dis = a s (Xs ~ Xi), (6) 

where a s is scale factor at source position. Xs, X\ are comoving 
distance from observer to source and lens. 

It is useful to introduce the tangential distortion compo- 
nent 7+ and the cross component 7 x which is rotated by 45° . 
The tangential shear components include all information of 



lens, if a lensing mass profile is axisymmetric on sky. These 
two components are defined with 7j and 72 as 

( T+ \ = ( -cos2i7 -sin2?7 \ / ^ \ 

\ 7x J \ -sin2>/ cos2>/ J \ y 2 J ' K > 

where ?/ is the angle between axis-a and 0. The sign of tan- 
gential shear indicate that background galaxy shape becomes 
tangentially deformed for positive and radially deformed for 
negative with respect to the lens center. 

Weak lensing signal is small so that it is difficult to mea- 
sure masses down to low mass objects. However, by stack- 
ing lensing signals from many samples, we can reduce er- 
rors which limit to obtain information on the objects (e.g . 
iMandelbaum et alj2006l;IOkabe et alj2010l;IOguri et al.l2012h . 
Therefore stacked lensing is a powerful tool for measuring 
average dark matter distributions of targets. The stacked 
lensing techniq ue is also used to constrain cosm ological pa- 
rameters (e.g. lOguri & Takadall201ll;lRozo et alJl201oi In the 
weak lensing limit, the avera ge tangential shear profile at 9 j 
is related to the convergence I Bartelmann & Schneiderl l200ll) 
as 

(7+) (00 =K(9<0 i )-(K)(9 i ), (8) 

where (■ ■ ■ } denotes average value in the circular annulus, 
the first term is the mean convergence within a circular aper- 
ture of radius 0j defined as 

ic(0 < 00 = i L a de ' K ( e ')' (9) 

and the second term is the mean convergence in the i-th ra- 
dial bin. The average cross shear component (7x}(0 < 0/) 
must be zero in weak lensing. 

In stacked lensing method, the following statistical un- 
certainty on tangential shear should be considered 

^ + = ^void + ^SS + <4ape' ( 10 ) 

where c" 2 oid is the statistical error coming from the deference 
in the structure of each void used in stacking analysis, c 2 g g is 
the error coming from large-scale structure along the line-of- 
sight and cr 2 ha p e is the shape noise error coming from the in- 
trinsic ellipticity of galaxies used for weak lensing measure- 
ments. Considering errors coming from structure of voids 
and LSS, the covariance matrix between i-th and j-th radial 
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bin is estimated as 



Hi 



Uv+M) ~ (7+(0i)»(? + ,;(0j) - <7+(0j))) 
j 

N v 



(11) 

where 7 + ;(0i) is the value of tangential shear in the 1-th 
void at i-th radial bin and N v is the number of voids. In 
addition to this, the error from intrinsic ellipticity shown in 
[Van Waerbeke et al. (2000) is determined as 



(12) 



sha P e N v n„S' 



where a e is the rms amplitude of the intrinsic ellipticity dis- 
tribution, n g is the number density of galaxies and S is the 
area of a bin. Throughout this paper we adopt a e = 0.4 
and tig = 30 arcmin~ 2 assuming survey parameters in the 
HSC survey. Shape noise becomes the dominant error at only 
small scale, on the other hand the other errors becomes dom- 
inant at large scale (sec !6.5t . We ignore other observational 
systematics such as imperfections in the telescope and elec- 
tric devices of detectors. 



The m ass is related to convergence iBartelma nn & Schneiderl 
2001) as 



m(< e) = J de' 2 K(e')\j(e' - e), 



(17) 



where V(9) is weight function and we adopt U(0) = 1 in 
our calculation. When we estimate void masses which are 
removed from underdensity regions, we calculate those only 
for the region with k{9) < 0. For the reconstruction from 
tangential shear, more careful calcula tion called ^-statistics 
jFahlman et al.lll994lClowe et alj200d) is needed. 

Furthermore the mass is related to the average mass 
density p and void radius R 



, 4 , 4 , 

M = - nW'p nW'p 

3 r 3 p 

= \\ntf{p-p)\ 



(18) 



= fR7 I 6 \, 
where we assume the mass density p is constant in the region. 



2.2 Signal to noise ratio of stacked lensing 

A useful way to quantify the observability is to calculate 
the total Signal-to-Noise (S /N) ratio over th e considered 
scale. Our calculation follows lOguri & Takadal j201ll) . The to- 
tal S/N is expressed as 

(|f ) 2 = E 7+ (ft) [cov( 7+ (9i), 7+ (0, ) )] -1 7+ (01 ), (B) 

where the summation indices i, j run over radial bins. For cal- 
culation with the tangential shear we estimate S/N for both 
situations that covariance matrix include and does not in- 
clude the shape noise. For calculation with convergence we 
can replace 7+ with k in eq[l3] The covariance matrix does 
not include the shape noise for this case. 

The S/N is related to the accura cy of the mass estimat e. 
More specifically, they are related as dOguri & Takadall201ll) 

„ ( £ V 1 (*n±y\ (i 4 ) 

Therefore S/N is proportional to the accuracy for constrain- 
ing the mass of the object. 

2.3 Mass reconstruction 

It is sometimes useful to derive a model indepen dent infor- 
mation of the mass from weak lensing data (e.g.|Okabe et alj 
l20ld: ISquires et all Il996l and ISquires et all Il997|) . Conver- 
gence in each circular annulus of radius 6 is related to the 
projected mass density 

m = as) 

*-*cr 

where Z cr is critical projected mass density defined as 

Ecr= 4^Gf^- (16) 



3 SIMULATION 

3.1 Raytracing simulation 

We use a l arge set of ray-tra cing simulation which are car- 
ried out in ISato et alj j2009l) . To carry out N-body simula- 
ti on, they use th e paralle l Tree-Particle Mesh code Gadget- 
2 <SpringellPi005] and see ISato et al]|2009l for simulation de- 
tails). The simulation employ 256 3 particles in a cubic box of 
240/z -1 Mpc on a side. Cosmological parameters are assumed 
as hubble parameter H = 73.2 km/s/Mpc, matter density 
n m = 0.238, baryon density n b = 0.042, dark energy density 
Q A = 0.762, equation of state parameter w = — 1, spectral 
index n s = 0.958 and variance of density fluctuation with 8 
h^ 1 Mpc ug = 0.76. T hese cosmological parameters follows 
WMAP 3-year results iSpergel et alj|2007h . 

By using results of N-body simulation, they carry out 
ray-tracing s imula tion with the algorithm dev eloped in 
IWhite & FIul teOOOl) and lHamana & Mellierl fcOOll) . They set 
source redshift z s at z s ~ 1.0 and make 200 realizations. The 
FOV in each run is 5 x 5 degree 2 . Therefore the total area is 
5000 degree 2 . 

Here we ou tline the ray-tracing simulation briefly, see 
ISato et all 12009) for more details. To obtain lensing signal, 
distance between observer and sources is divided into sev- 
eral regions to create set 19 lens planes separated by UOh^ 1 
Mpc in comoving distance. Then the three dimensional den- 
sity fluctuation in each each region is projected onto each lens 
plane. Next they trace 2048 2 rays backwards from the ob- 
server to the source plane. The initial ray directions are set on 
2048 2 grids with a grid space of 0.15 arcmin. Raytracing re- 
alizations of underlying density field are made by randomly 
shifting the simulation boxes assuming periodic boundary 
conditions. The lens planes from the same box are shifted in 
the same way in order to maintain the clustering pattern of 
mass distribution in one box. 
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3.2 VoidFinder 

For searching voids in the simulation data, we use 
the code called Void Finder which is publicly available 
([Foste r & Nelson 200 p). Th e vo id finding algorith m is based 
on lHoyle & Vogeleyl d2004h and lHoyle etaLM2005l) . 

At first galaxies in halo catalogue are separated into wall 
galaxies which are closer than the parameter R3 and other 
galaxies called field galaxies. The average distance to the 
third nearest galaxy D3 and the standard deviation 03 are cal- 
culated from halo catalogue, and R3 is defined as 



R 3 = D 3 + \cr 3 , 



(19) 



where A is dimensionless free parameter. A and minimum 
void radius £ affect the number and size of voids (see sec !6.U . 
In this paper, we ad opt A = 2.0, the same value used in 
iFoster & Nelson! i2009h . Wall galaxies are then set into cubic 
cell of R3 /2 on a side. The void finder finds the center of each 
empty cell. 

Secondly the maximal sphere are searched. Each empty 
cell is considered to be a part of void. The void finder sets 
it as a center of void, find the nearest galaxy and grow the 
sphere until a galaxy is on the edge. Next the void finder 
finds the vector connecting the center and the nearest galaxy, 
and mov away from the first galaxy along its vector and grow 
the radius of the sphere until the second nearest galaxy is on 
the edge. The vector connecting the center to the midpoint 
between the first and second galaxy is set. In a similar fashion 
to the second growth, the center along the vector is moved 
and the radius is grown until the third nearest galaxy is on 
the edge. In this algorithm voids are allowed to exist outside 
the boundary by extending it by 5 ~ 10 Mpc. 

Thirdly spheres which are found above are classified. 
The spheres are sorted from the largest to the smallest. The 
largest one is a void if the radius is lager than £. At first 
we consider only spheres whose radius exceed f . If a sphere 
overlaps other void by more than f> = 10% in volume, it is 
merged with that void. If not, it is defined as a new void. If a 
sphere overlaps more than two voids by more than ^3 = 2% 
in volume, it is rejected. Next we consider spheres whose ra- 
dius is smaller than f. Similaly if a small sphere overlaps 
more than two large spheres by more than ^3, it is rejected. If 
it overlaps one void by more than /3=50%, it is merged with 
the void. By this process, shape of voids may change into as- 
pherical shape. 

In this paper, we find voids that exist in the redshift 
range from z = 0toz = lin the simulation with this void 
finding algorithm. The void finder is applied to halos whose 
masses are more than 2.2 x 1O 12 M0 in the simulation. The 
total number of halos in the total FOV= 5000 degree 2 is 
1.2 x 10 7 within ^ z ^ 1. Therefore the number density 
is about 2000 degree - 2 . Then we use same pa rameters which 
are recommended in IFoster & Nelsonl 120091) except the pa- 
rameter of voids minimum radius in the process of void find- 
ing algorithm. We set minimum radius with 5 Mpc and 10 
Mpc and check the effect of changing the minimum radius. 

For stacked lensing analysis, we use voids with the radii 
between R — 15 ~ 45 Mpc and in the redshift range z = 
0.4 ~ 0.6 because lensing efficiency is expected to be high at 
the redshift range. The number of halos for selecting voids 
is 2.5 x 10 6 in the redshift range. Then for comparison, we 
also use voids with R = 15 ~ 40 Mpc in the redshift range 



z = 0.1 ~ 0.3. The number of halos used in this redshift 
range is 6.0 x 10 5 . 



4 PRESS-SCHECHTER THEORY 

For predicting the number count of voids, we use 
the modified Press-Sche chter (PS) formalism (see e.g. 
iKamionkowski et alj|2009h. 

In the PS theory (Press & Schechter 1974), regions which 
exceed the spherical collapse linear overdensity S c = 1.689 
become halos. The probability which has the linear density 
fluctuation S above S c becomes as 

1 (T-(R) ( &l 



P(5 > S C \M) 



-exp 



I2n S c " r V 2a 2 (R) 
where M is a halo mass defined with a radius R as, 



M 



47T 



-PR", 



(20) 



(21) 



(22) 



and a 2 (R) is the rms density fluctuation defined as, 

^(R) = / |^W 2 (/cR)P(fc), 

where W(fcR) is the Fourier transform of the window func- 
tion and P(k) is the matter power spectrum. Therefore the 
mass function n(M) is reduced as 



n(M)dM 

= 2p\P(S > S C \M) 
dlncr(R) 



•P(<5 > <S c |M + dM) 



2 j5_ 
nM 2 



dlnM 



«r(R) 



exp 



2cr2(R) 



dM. 



(23) 



In order to estimate the number count of voids based 
on the modified PS theory, we replace S c in eq|23]to linear 
underdensity \S V \. 

In our calculatio n, gaussian window function and BBKS 
iBard een et alj fl986) as transfer function are used. Void 
masses M v are estimated from the weak lensing signals with 
eq!18l For estimating the number count at a given redshift, 
the linear density fluctuation S v is evolved by using a linear 
growth factor. 



5 VOID MODEL 

Density profiles of voids are studied in prev ious works (e.g. 
ISheth & van de Weygaerll2004lPan et alj|2012h . These works 
indicate that in the under density region matter density is 
almost constant and there is a very sharp spike called a ridge 
near the edge of voids. 

In order to make a simple void model which include 
properties found in the previous works, we consider a spher- 
ically symmetric void model including ridges, named dou- 
ble top-hat model (fig©. One is for expressing voids, and 
the other is for expressing ridges. In this model, we set mass 
density p(r) as 



P(r) = { 



Pi 

Pi 





(r<Ri), 
(Rj < r < R 2 ) 

(r>Ri), 



(24) 



where pi and P2 are constant in each region. We assume that 
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the total mass between the void region and the ridge region 
should follow the mass conservation law. The mass between 
two regions is related as 

^nR^ip - Pl ) = i tt(R 2 3 - R! 3 )(p 2 - p). (25) 

From eq|25] p 2 is expressed with P \, Rj and R 2 

1 



P2 = P + ' 



■(Pl-P)- 



(26) 



1_(R 2 / Rl) 3< 

From eq|24]and eq|26] density contrast is reduced as 
r 6 (r<Ri), 
%) = ] FTO? (Ri<r<R 2 ), (27) 
I (r>R 2 ). 

An advantage of this double top-hat model is that its lensing 
properties can be computed analytically. By using eq|4]and 
eql27l convergence becomes as 



k(B) = 



3H -O m0 D ls D) /•D,v / i?^ 

D s a! 3 7-D, yV-fl 2 



2c 2 

+ 



dz{<50(Ri -r) 



1-(R 2 /R 3 )3 



e(r-Ri)0(R 2 -r) 



3H O m0 Di s Di 
2c 2 D sai 3 



2<5D 




02 



e\-e 2 



+ Je 2 - e 2 } 0(0! 



1 (02/d) 3 VM ^ ^ 

, ( &(e-e 1 )@(9 2 -e) 



e\-e 2 



i-(0 2 /0i) 3 



(28) 



where 0(r) is step function. The radii are related to angles 
r = Di^, (29) 



Dif?i, 



(30) 



where the induce i takes 1 or 2. 

In addition to the analytical expression of convergence, 
we can also calculate the tangential shear analytically. From 
eq(8] the tangential shear is reduce as 



(y+)(0) = R(<9)-(K)(e) 
2 

¥ jo 



[ de'e'Kte')- (k)(0). 

Jo 

From eq|28] the integration becomes 

/ de'e'Kte') = 

Jo 

' A(0 2 /0i) 3 (0?-0 2 )i -(0 2 - 2 )^ 



(31) 



3 l-(02/0i) 3 



0(0! - e) 



i)©( 



\ 3 l-(0 2 /0:) 3 
+0 ■ 0(0 - 2 ), (32) 

where A = 3H ^ 2 n m0 r^D? g Therefore with eq|28]and eq|32] 



P2 



Figure 2. Double top-hat void model. In r < Ri, p(r) equal to con- 
stant pi, and Ri < R < R2, p(r) equal to constant p^_. Relative am- 
plitudes between p\ and pi are determined from mass conservation 
law. Outside R2, p becomes zero. 



analytical expression of tangential shear is reduced as 

<?+>w = 

©(0! - e) 



30 2 1- (0 2 /0j) 3 
-{ie\ + e 2 )M-e 2 



'el - e 2 



a (2e 2 + e 2 )Je 2 -e 2 



J 30 2 1 - (0 2 /0i) 3 

+0 -0(0-02). 



0(0- 000(02 



(33) 



We assume the density contrast is equal to zero outside 
the second top-hat. Therefore the convergence and the tan- 
gential shear must be zero in that region. 

To extract parameters from simulations, we adopt the x 2 
method. The estimator x 2 is as follows 

X 2 = E (7+,model( e ;) - ?+,sim(0,)) 

x [cov(7 +;Sim (0,),7 +/Sim (0y))] 

) , (34) 

where ?+,model(fy) an d 7+,sim(^i) are tangential shear in the 
model and simulation at i-th radial bin. 



6 RESULT 

6.1 Number count of voids 

In order to check the dependence of the void number count 
on the parameter minimum radius £ in the void finder, we 
run void finder with £ = 5 Mpc and £ = 10 Mpc. 

Figure|3] shows histogram as a function of radius of 
voids for £" = 5 Mpc and £ = 10 Mpc for the redshift range 
^ z ^ 1. We can see the dependence on the minimum ra- 
dius £. The lager minimum radius is, the less number count 
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Figure 3. Number count of voids in the simulation for the redshift 
range < z ^ 1. x-axis is void radius, y-axis is number count of 
void in each bin. The dashed line and solid line histogram show the 
number count with minimum radius £ = 5 Mpc and £ = 10 Mpc. 
Error bar shows la with 200 realizations. 



of voids becomes. By using smaller minimum radius, we can 
obtain small voids which are not detected as voids with lager 
minimum radius. Therefore the number count of voids in- 
creases because the probabilit y of piling up with othe r voids 
increases. This trend matches Foster & Nelson] j2009|) . How- 
ever voids which have small radii in this calculation might 
not be true voids, because of the strong dependence with re- 
spect to the number count of voids on the value £. In the 
stacking analysis, we set the minimum radius £ = 5 Mpc. 
Figure|3] indicate that the number count of voids in the ra- 
dius range we consider in this paper (> 15 Mpc) does not 
depend strongly on the choice of the minimum radius. 



6.2 Direct estimate of void masses 

In order to estimate mass profiles and void masses, we stack 
lensing signals from voids selected by the void finder. From 
result of stacking and eq[17] we estimate the mass in each 
void sample. 

In sec !6.1l we obtained the number count of voids as a 
function of raius. We classify these voids into 6 bins based on 
the radii determined by the void finder. We select only voids 
that exist in the redshift range z = 0.4 ~ 0.6. We match cen- 
ters of voids determined by the void finder and stack those 
in each bin. The radii and number of voids used for stacking 
analysis are showed in table[TJ Points with errors in figured] 
and figure(5]show stacked convergence and tangential shear 
profiles in each radius. 

The void masses we estimate from the convergence pro- 
files and eq[l8] are presented in tablefT] The calculation is 
carried out only for the region k{9) < 0. We find that void 
masses that are directly estimated from the convergence pro- 
files are comparable to typical cluster masses (~ 10 14 M Q ).As 
we will see below, these void masses are smaller than true re- 
moved masses of the voids because of the effect of surround- 
ing ridges. 



Table 1. Details of stacking. Column (1): radius of voids determined 
by the void finder; Column(2): the number count of voids in the each 
radial bin and the redshift range zi = 0.4 ~ 0.6; Column(3): void 
mass calculated by integrating the k(9) < region. 

Radius [Mpc] Number of voids Mass from k(9) < [M G ] 



15- 
20 ■ 
25 ■ 
30 ■ 
35 ■ 
40 ■ 



20 

■■25 
30 

^35 
40 
45 



5246 
3892 
2446 
1400 
724 
320 



9.03 x 10 13 
4.45 x 10 14 
1.81 x 10 15 
3.06 x 10 15 
4.01 x 10 15 
3.26 x 10 15 



6.3 Fitting the convergence profile 

In order to compute convergence and tangential shear profile 
of voids analytically, we put a model named double top-hat 
model (sec|5). We consider the source at z s = 1.0, the lens 
at zj = 0.5. For our model, we set S, 9\ and 82 in eq|28l and 
eq(33]as free parameters to fit stacked lensing signals in our 
simulation. First we ignore shape noise and only consider 
noise from LSS and structures of voids in the process of fit- 
ting. Fitting is performed in the parameter range — 1 $J 8 ^ 0, 
0' sC 9 X < 100', 0' < 6 2 s$ 300'. We sum up x 2 over the ra- 
dius range ^ 9 ^ 90 arcmin for 15 ^ R ^ 30 Mpc, and 
< 9 sC 180 arcmin for 30 sC R < 45 Mpc. 

FigH] shows convergence profiles in simulation and the 
model. Model fitting is carried out for convergence profiles. 
We find that our simple model reconstructs lensing signals 
in simulation well. Best-fit parameters used in this plot are 
shown in table|2] Using these parameters, we estimate total 
void masses, i.e. the total mass removed from the void re- 
gion (table|2). Our model almost fits both convergence and 
tangential shear profile except the radius range 15 ^ R ^ 20 
Mpc. There are two reasons why fitting fails for the small- 
est void radius bi n. The first reason is a problem of finding 
voids. As we and lFoster & Nelson! J2009h showed, the num- 
ber of voids is sensitive to minimum radius £ in void finder. 
Therefore the void finder fails to find true voids which have 
small radii. The second reason is the violation of the mass 
conservation law. The mass conservation law is violated for 
simulated convergence profile presumably because of the 
large contaminations in the void catalogue at the small ra- 
dius bin. 

Compared with total void masses calculated from fitting 
of convergence profile with that from our model, masses es- 
timated from direct integration at k{9) < is a few times 
smaller. This is because observed convergence value is cu- 
mulative value from sources to an observer. Therefore when 
we calculate mass we should consider the effect from the 
ridge outside the void. Otherwise we underestimate total 
void mass. In addition, the size of underdense regions in the 
simulation and the model are smaller than radii defined in 
the void finder. This is probably because the halo catalogue 
we use for void finding does not exactly trace the underlying 
matter distribution. 
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Figure 4. Stacked convergence profiles estimated from simulation (points with errors) , the best-fit double top-hat model (solid line), assuming 
the survey area of 5000 deg 2 , the bin size AB = 6 arcmin and lens redshift zi = 0.5. Stacking analysis are carried out in each radius derived 
from the void finder; 15 < R < 20 Mpc (top-left); 20 ^ R ^ 25 Mpc (top-right); 25 < R < 30 Mpc (middle-left); 30 < R < 35 Mpc (middle-right); 
35 ^ R ^ 40 Mpc (bottom-left); 40 ^ R ^ 45 Mpc (bottom-rihgt). Best-fit parameters are summarized in table|2] 
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Table 2. Best fit parameters and derived total void mass. These parameters are determined from convergence profiles. Column (1): radius of 
void determined with the void finder; Column (2): density contrast; Column (3): radius of underdense region; Column (4): radius of overdense 
region; Column (5): void mass derived from the double top-hat model 



Radius in the void finder [Mpc] S [Mpc] [Mpc] Void mass [M, ;i ] 



15 - 


-20 


-0.390 


12.1 


25.5 


1.46 x 10 15 


20 - 


-25 


-0.194 


19.1 


27.5 


2.83 x 10 15 


25 r 


-30 


-0.244 


21.0 


63.0 


4.78 x 10 15 


30 - 


-35 


-0.424 


16.6 


109 


4.09 x 10 15 


35 r 


-40 


-0.316 


25.5 


59.8 


1.10 x 10 16 


40- 


-45 


-0.220 


27.3 


38.9 


2.73 x 10 16 



6.4 Fitting the tangential shear profile 

We also use tangential shear profile for finding best-fit pa- 
rameters. As we have done for convergence, we set same 
condition for searching best-fit parameters. 

Figure |5| shows tangential shear profiles in simulations 
and best-fit models. Table f5] shows best-fit parameters and 
total void masses determined from the fitting. Total void 
masses are estimated with eq[l8]as we have done in sec !6.3l 
We find that the best-fit parameters and the resulting total 
void masses are very similar to those obtained from fitting 
convergence profiles (table|2j, suggesting the consistency of 
our analysis. 



6.5 Signal-to-Noise ratio 

Using the stacked lensing covariance matrix we derived, we 
estimate S/N for convergence and tangential shear. For tan- 
gential shear, S/N is estimated in the case of with and with- 
out the shape noise. 

Fig[6]shows tangential shear profiles in both cases when 
the shape noise is added and not added to the covariance ma- 
trix. For illustrative purpose, points are shifted by 2 arcmin. 
Stacking analysis are carried out for voids in the radial range 
20 R ^ 25 Mpc and the redshift range z\ = 0.4 - 0.6. The 
shape noise is estimated from eq[12] The error from the shape 
noise is added to the diagonal elements of the covariance 
matrix. The effect from the shape noise becomes a dominant 
component of the error at small scale, because the number of 
galaxies in each radial annulus is proportional to the radius 
under the condition that radial bin size Ar is constant. Over 
the scale we consider, however, fig|6] shows that the shape 
noise is not a dominant component of the error. 

Total S/N is estimated for voids in the redshift ranges 
Z j = 0.4 ~ 0.6 and zj = 0.1 ~ 0.3 by using eqQU Tabled 
shows S/N in the redshift range zj = 0.4 ~ 0.6 for conver- 
gence, tangential shear without and with the shape noise as 
a function of radius. Table|5]shows S/N in the redshift range 
Zj = 0.1 ~ 0.3. TableH and tabled show that the staked 
lensing signals are detected at significant level (S/N^ 5). 
Adding the shape noise degrades S/N, but the effect is rela- 
tively minor. Compared with total S/N in the redshift range 
zj = 0.4 ~ 0.6, total S/N in the redshift range zj = 0.1 - 0.3 
are degraded. The error coming from structure of each void 
increases because the number of voids decrease as the red- 
shift range becomes low. However, S/N for individual void 
in the redshift range x\ = 0.1 ~ 0.3 is higher than that in the 
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Figure 6. Crosses and circles with errors show the tangential shear 
profiles for the cases which include the shape noise and not. The 
staking analysis are carrie out for voids in the radial range 20 ^ R $J 
25 Mpc. We assume a HSC-type survey, n g = 30 arcmin 2 , cr e = 0.4 
and FOV=5000 degree 2 . 



redshift range zj = 0.4 ~ 0.6 because of the larger apparent 
sizes of voids at lower redshifts. 



6.6 Comparison of the mass function of voids 

In order to predict the mass function of voids analytically, we 
determine the linear density fluctuation S v from comparison 
of mass functions between the simulation and the modified 
PS theory (sec[4}. In sec l6.3l and l6.4l masses of voids are es- 
timated for each radius. From this and sec !6.1| the number 
count of voids is obtained as a function of mass. We set S v as 
a free parameter. Fitting is performed in the parameter range 
0.1 sC |<5 V | s; 1.0. 

In order to obtain the mass function in each radius from 
the simulation, we use masses estimated from the tangential 
shear profiles (sec !6.4t . Because the inaccuracies of mass es- 
timation from weak lensing signals affect the mass function, 
we fit masses with the function as 



M = A R 



(35) 



where Aq is constant and set as a free parameter. Fig|7|shows 
the relation between masses and radius estimated from void 
finder and the best-fit curve. From the least-square method, 
Aq is determined as 



1.32 x 10 



11 



(36) 
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Figure 5. Tangential shear profiles estimated from simulation (points with errors), the best-fit double top-hat model (solid lines) in each radius. 
Same assumptions in fig|4]are adopted. Stacking analysis are carried out in each radius derived from the void finder; 15 ^ R ^ 20 Mpc (top-left); 
20 < R < 25 Mpc (top-right); 25 ^ R ^ 30 Mpc (middle-left); 30 ^ R ^ 35 Mpc (middle-right); 35 ^ R < 40 Mpc (bottom-left); 40 < R < 45 Mpc 
(bottom-rihgt). Best-fit parameters are summarized in table|3] 
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Table 3. Best-fit parameters and derived total void mass. These parameters are determined from tangential shear profiles. Column (1): radius of 
void determined with the void finder; Column (2): density contrast; Column (3): radius of underdense region; Column (4): radius of overdense 
region; Column (5): void mass derived from the double top-hat model 



Radius in the void finder [Mpc] S [Mpc] [Mpc] Void mass [M, ;i ] 



15 - 


-20 


-0.398 


12.2 


25.3 


1.53 x 10 15 


20 - 


-25 


-0.170 


19.1 


27.6 


2.50 x 10 15 


25 r 


-30 


-0.250 


19.0 


109 


3.63 x 10 15 


30 - 


-35 


-0.403 


16.8 


40.0 


4.00 x 10 15 


35 r 


-40 


-0.290 


25.5 


59.8 


1.01 x 10 16 


40- 


-45 


-0.183 


27.3 


39.5 


7.80 x 10 15 



Table 4. Signal-to-noise ratio, integrated over angular scales considered in the redshift range z = 0.4 ~ 0.6. (1): radius of void determined with 
the void finder; Column (2): the number count of voids in the each radial bin and the redshift range z = 0.4 ~ 0.6; Column (3): S/N derived 
from convergence; Column (4): S/N derived from tangential shear without the shape noise; Column (5): S/N derived from tangential shear 
with the shape noise; Column (6): S/N for one void. It is estimated from information of tangential shear with shape noise and the number of 
voids (table[lj 



Radius [Mpc] Number of voids 



S/N 
(convergence) 



S/N (tangential shear 
w/o shape noise) 



S/N (tangential shear 
w/ shape noise) 



(S/N)/VN—d 
(w/ shape noise) 



15 ■ 
20 ■ 
25 . 
30 . 
35 ■ 
40 - 



■ 20 

■ 25 
- 30 
■35 

40 
45 



5246 
3892 
2446 
1400 
724 
320 



11.7 
9.60 
8.65 
8.17 
8.08 
7.63 



11.0 
8.25 
7.55 
8.19 
7.57 
7.89 



9.18 
6.97 
6.51 
6.76 
6.25 
5.39 



0.127 
0.112 
0.132 
0.181 
0.232 
0.301 




20 25 30 
radius [Mpc] 



35 40 45 



Figure 7. Mass-radius relation, x-axis is radius obtained from void 
finder, y-axis is mass estimated from our model with tangential shear 
profiles. Crosses show masses obtained from weak lensing signals 
and the dashed line plots best-fit curve with eql35l 



Using eq|35] we divide the number count obtained from sim- 
ulation by masses included in each radius range. The number 
count of voids is used in the redshift range 0.4 < z ^ 0.6. 

Figure H| shows the mass function in the simulation and 
the best-fit curve of the mass function in the modified PS the- 
ory. Mass functions of voids in simulation are reproduced 
with S v = —0.35. For checking the consistency, we also do 
same works for the case of convergence. The best-fit linear 
density becomes as S v = —0.5 in that case. These values 
are smaller than the value predicted from the spherical top- 
hat model (Sheth & van de Weygaer3l2004]) . From our result, 
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Figure 8. Mass function of voids in simulation and the modified PS 
theory. Points with errors show mass function of voids derived from 
the simulation. Masses are estimated from tangential shear profiles. 
The number count of voids in the redshift range 0.4 ^ z ^ 0.6 are 
used. Error bar shows la of the number count with 200 realizations. 
The dashed line shows the best-fit curve in the modified PS theory. 
The linear density fluctuation S v = —0.35 is adopted. 



however, the linear density parameter is sensitive to the mass 
estimation reconstructed from weak lensing signals. For bet- 
ter understanding, we should study structures of voids and 
ridges more carefully. 
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Table 5. Signal-to-noise ratio, integrated over angular scales considered in the redshift range z = 0.1 ~ 0.3. Column (1): radius of void deter- 
mined with the void finder; Column (2): the number count of voids in the each radial bin and the redshift range z = 0.1 ~ 0.3; Column (3): S/N 
derived from convergence; Column (4): S/N derived from tangential shear without the shape noise; Column (5): S/N derived from tangential 
shear with the shape noise; Column (6): S/N for one void 



Radius [Mpc] Number of voids 



S/N 
(convergence) 



S/N (tangential shear 
w/o shape noise) 



S/N (tangential shear 
w/ shape noise) 



(S/N)/VN— d 
(w/ shape noise) 



15 ■ 
20 . 
25 ■ 
30 ■ 
35 . 



■ 20 
< 25 
> 30 
■35 
40 



1798 
1097 
523 
252 



6.49 
5.73 
7.36 
5.65 
6.93 



5.69 
7.21 
5.09 
6.62 



5.22 
6.70 
4.43 
5.67 



0.144 
0.158 
0.293 
0.279 
0.605 



7 CONCLUSION 

In this paper, we have presented the observability of voids 
with stacked weak lensing in next generation surveys such 
as HSC. We select voids with a void finder from halo cata- 
logue made from a large set of N-body simulations. We have 
stacked convergence and tangential shear profile for them 
with the full ray-tracing simulations to obtain an average 
dark matter profile (black in fig|4|and fig|5). From our stack- 
ing analysis, we have seen both structures of void and ridge 
outside the void. At small radii, the void finder cannot find 
true voids if the minimum radius £ is set to small value. This 
resul t is consistent with the previous paper I Fos ter & Nelsonl 
2009). In oder to predict the mass function of voids analyti- 
cally, we compare the mass function of voids in the simula- 
tion with the modified PS theory and found that the mod- 
ified PS theory can predict the mass function of voids ob- 
tained from the analysis for tangential shear profiles when 
the threshold linear density fluctuation is set to S v = —0.35. 
In addition, we have estimated total S/N, which indicates 
that stacked lensing signals from voids can be detected at sig- 
nificant level (S/N^ 5) for the 5000 degree 2 area, even if the 
error coming from the shape noise is added. The shape noise 
error becomes a dominant component at only small scale. 

To fit lensing signals in stacked analysis, we have con- 
sidered a simple void model called double top-hat model. 
Our model fits both profiles of convergence and tangential 
shear in simulation very well (fig|4|and figO- Estimated to- 
tal void masses from this model were M = 10 14 ~ 10 16 Mrv,, 
which were a few times larger than masses derived from di- 
rect integration of the convergence profile at k{6) < 0. We 
have confirmed that the dense ridges outside voids affect 
profiles of weak lensing signals and the mass estimation of 
voids. Therefore we have to take account of this effect for in- 
terpreting stacked weak lensing signals around voids. The 
underdense region is smaller than radius derived from the 
void finder, presumably because the halo catalogue used for 
finding voids does not trace true dark matter distributions 
directly. 

Our work has demonstrated the observability of voids 
with the stacked weak lensing method. We have also 
suggested a simple model to interpret the lensing signals. 
In this work, the average number density of halos is 200 
/degree 2 /Az = 0.1 in the redshift range z = 0.4 ~ 0.5. In 
the SDSS-HI Baryon Oscillation Spectroscopic Survey (BOSS) 
{Dawson e t al. 2012), the average number density of lumi- 
nous red galaxies (LRG) is about 50 /degree 2 /Az = 0.1 in 



the same redshift range. From eq[l9] void radii we can find 
mainly depend on the average distance between galaxies. 
In the simulation, it becomes about 6.86 Mpc. In BOSS, it 
becomes about 10.9 Mpc, which is about a factor of two 
larger than the average distance used in this work. However, 
it would be possible to find voids from the BOSS data which 
have large radii, and adopt techniques described in this 
paper using weak lensing measurements in e.g. HSC data. 
These techniques can al so be applied to n ext generation 
surveys such as Eu clid jLaureiis et alJl201l|). Dark Energy 
Surv ey (DES) {The Dark Energy Survey Collaboration! 
2005) and Large Synoptic Survey Telescope (LSST) 
JLSST Science Collaboration et alj2009n . 

W hen we were writ ing this paper, we came across a pa- 
per by lKrause et alJ te012l) which also studied stacked weak 
lensing signals around voids. The main differenc e of their 
paper from our analysis is that lKrause et alJ J2012I) assumed 
analytic mass profiles to estimate S/N, while in this paper 
we present realistic stacked lensing profiles and the error co- 
variance with ray-t racing of N-body simulations. In addition, 
iRrause et al.l j2012l) considered voids with smaller radii of 
^ 15fe _1 Mpc for which our void finder does not work well. 
Nevertheless, by extrapolating both the results we argue that 
they are broadly consistent with each other. 
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